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TECHNICAL NOTE 3712 


BOUNDARY LAYER BEHIND SHOCK OR THIN EXPANSION 
WAVE MOVING INTO STATIONARY FLUID 
By Harold Mirels 


SUMMARY 

The boundary layer behind a shock or thin expansion wave advancing 
into a stationary fluid has been determined. Laminar and turbulent bound- 
ary layers were considered. The wall surface temperature behind the wave 
was also investigated. The assumption of a thin expansion wave is valid 
for weak expansions but becomes progressively less accurate for strong 
expansion waves. 

The laminar-boundary-layer problem was solved by numerical integra- 
tion except for the weak wave case, which can be solved analytically. 
Integral (Kerman -Pohlhaus en type) solutions were also obtained to provide 
a guide for determining expressions which accurately represent the numer- 
ical data. Analytical expressions for various boundary-layer parameters 
are presented which agree with the numerical integrations within 1 percent. 

The turbulent -boundary-layer problem was solved using integral meth- 
ods similar to those employed for the solution of turbulent compressible 
flow over a semi-infinite flat plate. The fluid velocity, relative to 
the wall, was assumed to have a seventh-power profile. The Blasius equa- 
tion, relating turbulent skin friction and boundary-layer thickness, was 
utilized in a form which accounted for carpreBsibility. 

C ons ideration of the heat transfer to the wall permitted the wall 
surface temperature, behind the wave, to be determined. The wall thick- 
ness was assumed to be greater than the wall thermal-boundary-layer 
thickness. It was found that the wall temperature was uniform (as a 
function of distance behind the wave) for the laminar-boundary-layer case 
but varied with distance for the turbulent -boundary-layer case. 


INTROIUCTION 

If a shock or exp ans ion wave advances into a stati o n ar y fluid bounded 
by a wall, a boundary -layer flow is established along the wall behind the 
wave. boundary-layer flow is often important in studies of phe no mena 



2 


NACA IN 3712 


involving nonstationary waves (e.g., shock tube studies, initiation of 
detonation studies, etc.)* * 

In references 1 to 3, the laminar boundary layer behind a shock 
wave is studied. The purpose of this paper is to extend reference 1 to 
include the case of thin expansion waves (i.e ., expansion waves of zero 
thickness) and the case of turbulent boundary layers. This extension 
was motivated primarily by a desire to provide boundary-layer informa- 
tion for studies of shock attenuation in a shock tube (e.g., ref. 4). to 

co 

• The assumption of a zero-thickness expansion wave (referred to as 
a "negative shock" in ref. 4) becomes progressively less accurate as the 
strength (pressure ratio) of the expansion wave increases. Thus, the 
applicability of the boundary-layer results in shock tubes wherein strong 
expansions occur may be questioned. Nevertheless, the assumption of a 
zero-thickness expansion wave will be retained for ranges of the param- 
eter Uy/ug (symbols defined in appendix A) corresponding to strong ex- 
pansion waves. The reason for retaining the assumption is twofold. 

First, a qualitative estimate of the boundary layer in the conBtant- 
pressure region behind strong expansion waves is obtained. It may be , v 

possible to modify these results to give better agreement with more ac- 
curate estimates of the boundary layer behind a strong expansion wave. 

Second, the results are directly applicable to other physical problems 
such as (l) the boundary layer on a semi-infinite flat plate with tan- 
gential blowing and (2) the boundary layer behind a planar flame front 
propagating in a stationary combustible mixture. 

An independent investigation of the boundary layer behind an ex- 
pansion wave, talking into account the finite thickness of the expansion 
wave, is currently in progress at the MCA Langley laboratory. 


ANALYSIS 

A shock or expansion wave is considered to move with constant veloc- 
ity, parallel to a wall, into a stationary fluid (fig. l). The boundary 
layer along the wall behind the wave is to be determined. In the case 
of an expansion wave, finite thickness of the wave will be ignored. Both 
i rnninn.T* and turbulent boundary layers will be considered. 


Coordinate Systems 

Let x,y be a coor dinat e systemjfixed with respect to the wall and 
let u,v be velocities parallel to x,y as indicated in figure l(a) . 
The flow is unsteady in the x,y-coordinate system. Let x,y represent 
a coor dinat e system moving with the wave (fig. l(b)). The velocities 
parallel to the x- and y-coordinates are denoted by u and v, respec- 
tively. In this coordinate system, the flow is steady. 
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_Assume that at time t = 0 the two coordinate systems coincide. 

If is the velocity of the wave relative to the wall, then x and 

x are related by * = J Iba axial velocities are related by 

u a u - u^. Bote that the wall moves with velocity u^. « -u^ in the 

x,y-coordinate system. The transformation relating the two systems' can 
also be expressed as 

X « X - u w t| 


u 3 u - u w 
V a V J 

Equations defining conditions across the shock and expansion waves are 
summarized in appendix B. 



Inrnimr Boundary Layer 

The Prandtl boundary-layer equations apply for the flow in the vi- 
cinity of the wall, except at the base of the wave where the boundary- 
layer assumptions break down. By assuming laminar flow and dp/dx « 0, 
the boundary-layer equations for x > 0 are 


^u + c 
ax c 


(2a) 

du . chi 

1 h f ba\ 

(2b) 

U dx Sy * 

* p $y v / 

/ St . &r\ 


(2c) 

P 83 

pBT 

(2d) 

The boundary conditions for x > 0 

are 


u(x,0) => u^ 
v(x,0) = 0 

u(x,») « u e ^ 

► (3) 

E(x,0) a 

T(x,«) a T e 
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These equations can he transformed to a form suitable for numerical in- 
tegration. In accordance with, reference 5, it will he assumed that |i 
and k are proportional to T and that c^ and cr are independent 
of T. These state properties will he defined to have the correct numer- 
ical values at the wall. The use of a mean reference temperature, rather 
than a wall reference temperature, is discussed in appendix C. 


From, equation (2a), a stream function i|r exists so that 


p w 

£1 ■ 



(«■) 


Define a similarity parameter 



and assume that i|r can he expressed in the form 

+ - ( 6 ) 

The velocities are then 



Assume that p. and k are proportional to T and choose the constants 
of proportionality so that p and k have the correct numerical values 
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at the wall; similarly, assume that c and a axe constant at their 
wall values. Thus, ^ 

n -sr I 

w 


^w 

w 


C S3 C 

P P^W 


w 


J 


The momentum equation can now be written 

f »" + ff" » 0 
f(0) a 0 

f'(o) a uyu e 

f '(«.) a 1 J 

and the energy equation becomes 


(t j * V® 




a u‘ 


w e 




(f") 


n\2 


T(x,0) ^w 




TCx,”) 


• T 


T 


a 1 


J 


( 8 ) 


(9) 


(10) 


Since equations (10) are linear, T can be expressed as the linear super- 
position of the solution for zero heat transfer plus the effect of heat 
transfer. That is, 


€- i+ £- x ) : 


A 2 + _*£!.(,) 

®e«p,» V T e .%! 


(H) 
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■where 


^r . J* w \ 2 u e r (Q) 

^e^' \ u e ) 

and r(rj) and s(tj) are defined, respectively, by 


2 0 

r" + cx w fr« = — - V o (f ") 2 


(H 

r(«) = r 1 (o) « 0 


and 


( 12 ) 


(13) 


s" + C^f s ' a 0 

s(O) a 1 
s(«) a 0 

Equations (13) and (14) can also be expressed in quadrature form: 


(14) 
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(16c) 


(16a) 


where q equals the heat transfer in the +y-direction . Equati on (l6a) 
results from the expression Vg/u e = d5*/dx., which can he deduced from 
the continuity equation and the definition of 5*. 


Equations (9) to (l5) can he integrated analytically for 
l^y/ug - l| « 1 (i.e. j weak waves). The integrations were conducted 
in reference 1 and the results are summarized as follows: 
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f o rd ""^ + 0 S‘ 1 ) 



Cl7f) 


(I7g) 

(I7h) 


(171) 


Equations (15) are integrable for a w =* 1. 


The results show 



(18a) 


(18b) 


(18c) 


(l8d) 


Equations (9) have been integrated numerically (using an IBM card- 
programmed electronic calculator) for u^/ug » 1.5,- 2, 3, 4, 5, and 6, 
and the results are presented in reference 1. The range 1< u„./u e « 6 
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covers the entire range from, very weak to very strong shock waves. 
Equations (15) are also evaluated in reference 1 for a v - 0.72 and 
u^/rig a 2, 4, and 6. In order to extend the numerical calculations to 


the case of a thin expansion wave, equations _ (9) have teen integrated for 
u^/ug =» 0, 0.25, 0.50, and 0.75 (table I) . The case u^/ug =* 0 is the 

well-known Blasius problem for flow past a semi-infinite flat plate and 
is included herein for the sake of completeness. Equations (15) were 
integrated for Uy/ug s» o, 0.25, 0.50, 0.75 and cr w sa 0.72 (table l). 

Prom appendix B, the limiting value of u^/ug for strong expansion waves 


, ”w 
is — 
u„ 


r + i* 


These calculations were made using the integration tech- 


nique described in appendix B of reference 6. The results of reference 1 


for u^/ug =» 2, 4, and 6 are also included in table I. 


Integral solution . - The numerical integrations of equations (9) 
and (15) have the limitation that only specific values of i^/ug and 
cr w can be treated. The analytic dependence of the boundary-layer pa- 
rameters on u^/ug and a w are not explicitly defined by these specific 

solutions. In order to obtain approximate analytic solutions, an inte- 
gral (Karman-Pohl hausen) method of approach will now be used. The re- 
sults will be used later in this section to obtain analytic expressions 
which accurately represent the numerical data presented in table I of 
this report. 


Integrating equations (9) with reBpect to q, in the range 
< T) < », gives the integral form of the momentum equation 


f"(0) 


■s: 


f’(l - f 'Jdq 


(19) 


Let y =3 8 be the edge of the boundary layer from a Pdhlhausen point of 
view. The corresponding value of q is desig n ated Tig. Define 


£ = ri/rjg . The velocity f 1 
nomial in £ in the range 
tions on f 1 noted in equations (9) plus the boundary conditions 
(f"V S3 o, in order to satisfy equations (9), 

b°(j 


is now represented by a fourth-degree poly- 
£ < l. 1 Using the two boundary condi- 


and 


^References 3 and 7 also used integral approaches employing fourth- 
degree polynomials in studies of the laminar boundary layer behind moving 
shocks. In reference 8, a normalized error function was used instead of 
the fourth-degree polyn omial - This resulted in very good agreement with 
the exact numerical integration of the corresponding boundary-layer 
equations . 
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(f") *» (f = 0, in order to get a smooth joining with the ex- 

S” 1 5=1 

temal flow, gives 


for 0 < £ <. 1, 


for 1 < £, 




f * a 1 


Equation (l9) can he written as 


J 


; " (0) “ *jf 


f*(l - f)<l£ 


Substituting equations (20) into equation (2l), noting 


gives 


% 


2 


315 


i Utr 

1 74 + 115 ^ 
u e 


(20) 


( 21 ) 


( 22 ) 


uhich, together with equations (20), defines f • . The following quan- 
tities are of interest: 


f"(0) 


(‘■-S' 


74 + 115 — 
Ue 


315 


(23a) 


lim (q - f) 
rr** 



(1 - f«)dTl 


_ 6 _ 

10 



(23b) 


Equations (23) are compared, in table XX, with the exact integration of 
equations (9) for u = 0 , 1, and 6. It is seen that equation (23a) 
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agrees within 4 percent with the exact integration. Equation (23b) 
agrees within 2 percent at u^/ug » 0 and within 10 percent at 

= 6. Table H serves to indicate the accuracy of the integral 
method of solution of the momentum equation. 


Substitution of equations (20) and (22) into equations (15) gives 

r 3 &( o v ,t) 


where 


/ r dii - T]glCa w ) 

p 00 

/ S dl) « Tlgj(o w ) 

'-'O 


8 1 (o) =* K(a w )/rig 
*(£} s (1 - £)*( 2£ + 1) 



(24a) 

(24b) 

(24c) 


(24d) 

(24e) 

(25a) 


<*C<V£) = C25b) 


(25c) 


(25d) 



(25e) 


(25f) 
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Equations (24) define the functional dependence, on cr w ,£ and u^/ug, 

of the various quantities noted to within the accuracy of the integral 
method. Equat ions (25) can he integrated analytically for ■ ff *» 1. For 
other Prandtl numbers, numerical integrations, or approximate analytical 
integrations would be required. 


Equations (24c) to (24e) indicate that the functional dependence of 

n« r\ co 

/ r drj, / s dq, and s'(o) on is independent of u^/ug. 

Comparison of equation (24d) with the exact solution for s’(0) at 
u^/u e = 1 (i.e., eq. (l7h)) Indicates that K(a w ) should be proportional 

to cr^ 1 / 2 . However, comparison of equation (24d) with the exact solution 
tion for u^/u e = 0 (ref. 5) indicates that K(a w ) should be proportional 


to cr 1 / 3 for a in the range 0.6 < cr T , < 1. The discrepancy is due to 
w w 

to limitations of the Integral method. Hence equations (24) and (25) 
would have to be modified if greater accuracy is desired. Possible mod- 
ifications are indicated in the next section. Comparison of equations 
(24c) and (24d) with equations (l7f) and (l7i) indicates that, for u„/u e 

near 1, l(a w ) is a constant while J( G W ) varies directly as l/ 


interpolation formulas - - Equations (23) and (24) can be used as a 
guide to obtain analytic expressions which accurately represent the data 
obtained from the exact numerical Integrations of equations (9) and (15). 
Thft various boundary-layer parameters are assumed to depend on u^/u e 
and a._ an follows : 

W 



(26a) 


(26b) 




D 5 +D 6K/ U e) 


(26c) 
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s di] 


i 


1 + Co — 
8 ^ 


(a ) 
v w' 


^eCSr/ue) 


(26d) 


s«(0) = C c 


11 + C 10 u e ^ 


u7, ,V D 10^V u e) 


(26e) 


r(o) - (<* w ) 


D ii +D i2K/ u e ) 


(26f) 


The case of an expansion wave and a shock wave are distinguished by 
writing separate fc 
1< u^/Ug-^ 6. The 

for 0 < u^Ug 4 1, 


writing separate formulas for the ranges 0 <j u^/ug^ 1 and 
K u^/Ug^ 6. The results are, 


— ^ - - 0.470 ill + 1.887 — 

1 


13m ILS-1 ... , -..1-2 1? 


-s i 


1 + 1.326 — 
u e 



1.686 


Ml + 1.234 — 

I 




0.22(l-u w /u e ) 


1.217 




1 + 1.326 




(o ) 


-0.36-0.14 Cu^/Ug ) 


u„ 


s' Co) 


0 . 35 + 0 . 15 ( 11 ^/^) 


-0.470 /\jl + 1.887 ^ (a w ) 

0 .50-0 .13(13^/^) 


r(0) - (o w ) 


(27a) 


(27b) 


(27e) 


(27d) 


(27e) 

(27f) 
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for 1 < £%/%) < 6, 



» 0.489 A 1 + 1.665 — 

1 _ \ u e 

u„ * 


H — f 1.155 


11 m 

T]-X» i _ ~w 


1 ^ 

1 u~ il/l + 1-022 — 

U e ^ u e 


r dtj 


1.569 


il/l + 0.933 ^2£ 

If 




0.045(1-11^/113) 



dt] 


1.134 


/l| 1 + 1.022 - 

If 


% 

. 

e 


K) 


-0.47-0.029(u„/u } 


1 


-0.489 A 1 + 1.664 ^ (o w ) 


Us. 0.48K).022(u w /u e ) 


r(0) 


0.39-0.023 (Uy/ u e ) 


(27a'} 
(27b ' ) 
(27c ' ) 

(27d») 

(27e') 

(27f') 


Equations (27) were obtained by matching equations (26), at the end 
points of the intervals 0 < (t^/ug) < 1 and 1< (u^/Ug) < 6, with 
the exact integrations of equations (9) and (l5) (table i) . Equations 
(27) agree with table X to within 1 percent. The accuracy of equations 
(27) can be estimated from, table HI, which compares equations ( 27 ) 
with table I for the intermediate points u^/ug =» 0.25, 0.50, 0.75, 2, 
and 4. 


Equations (27) are based on numerical solutions corresponding to 
a w a 1 and a w =» 0.72. They should give reasonable estimates for other 

values of o w of this order (i.e., 0.6 < l). An exact appraisal 

of equations (27) for a w other. than 1 and 0.72 requires further inves- 
tigation of the exact integrations of equations (9) and (15). Such an 
Investigation is considered beyond the scope of this report. 
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W«.t 1 surface temperature ♦ - Until now it has been assumed that the 
wall surface temperature is constant behind the shock or expansion wave. 
This can be shown to be true for a laminar boundary layer. The magni- 
tude of the wall temperature will now be determined, assuming that the 
wall thickness is greater than the wall thermal boundary layer. 2 


In the steady x, y-coordinat e system, the wall has a uniform veloc- 
ity Uy.. However, there is a thermal boundary layer in the wall which 
is s imilar to that in the adjacent fluid (fig. 2) . Assume that before 
the appearance of the wave, the fluid and the wall were in thermal, equi- 
librium at the temperature T^. 

When a q uant ity is discontinuous across y 53 0, the subscripts w 
and 0“ will be used to indicate quantities evaluated at y » 0 + and 
y « 0", respectively. Neglect the variation of state properties (such 
as k, p, etc.) in the wall and write the diffusivity of the wall as 

. By using boundary -layer assumptions, the energy equation 

0 " 

for the wall becomes 

» . V aV (28a) 

S ”v 8/ 1 



with the boundary conditions, for x > 0, 

T(x,0) « 0^. 

T(x, -®) a 0^ 


(28b) 


Equation (28a) can be obtained from equation (2c) by assuming the wall 
to he equivalent to a fluid with, p, « v » 0. Assume (to he verified 
later In this section) that T^ is constant. The solution of equations 
(28) is then 


*w - T 

*w- T b 



(29) 


^The problem of the w all surface temperature associated with the 
iBrninwr boundary layer behind a moving shock was solved, independently, 
in reference 8. The development in reference 8 is somewhat more de- 
tailed t han that given herein. The results are essentially in agree- 
ment with this report. 
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The heat transfer in the +y-direction at y = 0“ is 



( T b - T vV 

=3 ■ ■ — 

From equation (l6b), the heat transferred into the fluid is 
Equations (30 ) and (31) may he equated, yielding 

*w - * b . [C*r/*b) - l] 

% “ 1 + A 


(30) 


(31) 


=* A 



+ 0(A 2 ) 


(32a) 


where 


A 


-s 1 (o) . Fe a I°0- \ 

'yjt w 


(32b) 


Equation (32a) is solved for T by iteration (since A and T are 
functions of T w ) . The choice T w = T^, to start the iteration, is a 
good first approximation. Note that is a constant, which is con- 

sistent with and, therefore, substantiates the original assumption to 
this' effect . 


The ratio T r /T-^ can be expressed as a function of v^/ug and 

r (o) by use of equation (12) and equations (B3), (B4), (B7), and (B8) 
of appendix B. The results are (assuming c^ e = c ^ 
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for V u e > ^ 


5 

Ti. 


-! . (i . A 

\e ) / r + 1 _ % 


r (o) 


/ r + 1 _ SA 

u e V T " 1 u e/ 


(33a) 


%• 

= 2 ^S— » | (l^ - l) for y = 1.4; r(o) » 1 


6 - 

u e 


for uJu Q < 1, 

*- \ 
2 [ Q C u w/^ e) " 3](l 


- 1 (33b) 


* s K/%y‘ 


w - 2 



1 

d 

1 

1/7“ 

/*n \ 

2 

L 3 l^ - 2 J 

I *6 1 

t " W 


> for r = 1.4; r(o) = 1 


“V 


It should be noted that for T r /T b > 1 the heat transfer is from 
the fluid to the wall, whereas for < 1 the heat transfer is from 

the wall to the fluid. Thus, for the flow behind a shock wave (eg. (33a)) 
the heat transfer is always from the fluid to the wall. However, for the 
flow behind an expansion wave (eq. (33b)) the heat transfer is from the 
wall to the fluid for relatively weak expansion waves and from the f luid 
to the wall for strong expansion waves. When y =« 1.4 and r(o) « 1, 

the change in the direction of heat transfer occurs at P e /]?b = (2/3) 7 . 

This change in heat-transfer direction is a consequence of the fact that 
the stagnation temperature (relative to the wall) behind a thin p- giRna -inn 

2 x_ 

wave is less than T^ for Pg/Pt, 

2 L. 




and is greater than T v 


for Pg/Ph 


(mr- 


The boundary-layer recovery temperature T t 
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equals the free-stream stagnation temperature (relative to the wall) 
for r(o) => 1 and is somewhat less than the stagnation temperature for 
r(0) <1. If the finite thickness of the wave is taken into account 
for a very strong expansion wave, it is expected that the heat transfer 
will he from the wall to the fluid for the early stages of the expansion 
and from the fluid to the wall for the final stages of the expansion. 


The order of magnitude of (T^ - 

Equations (l7h), (27e), and (27e') indicate 
Therefore, A can he approximated as 


will now he determined. 


-s’( 0 ) 




u e 

— si. 
"w 


A 


fit 



(34) 


If k is proportional to T and Cp Is constant, then the value of 
k/-^/a is independent of temperature for a given gas. Assume the fluid 
to he air and the wall to he cast iron (carbon, approx. 4 percent) hoth 
initially at a temperature 1^ » 530° R. Using the property values of 

reference 9 evaluated at 530° R gives A s 0.0004. Substituting equa- 
tions (33) into equation (32a) yields, with A =» 0.0004, r(o) * 1, and 
T =* 1.4, 

for u^/u e > 1, 

T ^ - s 0.0004 (35a) 



(35b) 


The wall will remain within 0.1 percent of its original temperature for 
shocks up to s 3 and for an expansion waves. The wall will remain 

within 1 percent of its original value for shocks up to s 9 and 
within 10 percent of its original value for Mp = 27. Similar results 

would he obt ain ed if other gases a.nri wall metals are used. The impor- 
tant parameter k w /kQ_ is small because of the relatively low heat con- 
ductivity of gases c omp ared with metals . The fact that the wall remains 
essentially at its original temperature (except possibly for extremely 
strong shock waves) justifies the disregard of the variation of wall 
state properties with temperatures in equation (28a) . 
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The thickness of the wall thermal, boundary layer Ap can readily 
be found. Define the edge of the thermal boundary layer to correspond 
with (l^ - Tj/tl^ - Ijj) = 0.99. Then, from equation (29), 

Ap - 3.6^xctQ_/u w . (36) 

It is of interest to compare Ap with the thermal -boundary-layer thick- 
ness in the fluid Stp. For convenience, consider the case of weak waves. 

Equations (ll) and (l7g), together with the previously defined criterion 
for the edge of a thermal boundary layer, give 



which is of the order of 1 for some air and metal wall combinations . 

The solution for the wall surface temperature presented herein is 
valid provided the wall thickness is (at least) greater than Ap. 


Turbulent Boundary Layer 

Previously, it was assumed that the boundary layer behind the wave 
was laminar. Now the turbulent -boundary-layer case will be treated 
using integral methods . The solutions will be obtained by extend in g 
empirical, semi- infinite flat -plate, boundary-layer data to the case 
where the wall is moving. 

Boimflm-y-iayer solution . - The integral form of the momen t um equa- 
tion is 



d6> 

dx 


(39) 


■where u/u represents the average velocity. As before, denote the 
boundary -layer thickness by 5 and define a similarity parameter 
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^ b y/s. It is necessary to express u/ug as a function of In 

studies of the turbulent boundary layer on a semi-infinite flat plate, 
it is customary to assume u/u g = ^V 7 for 0 < ^ 1. The use of the 

seventh-power law appears to agree with available data, even for high 
Mach number flows, and to result in reasonably accurate relations for 
skin friction (e.g., refs. 10 and 11 ) . In order to extend this profile 
to the case of a moving wall, it is assumed that the fluid velocities 
relative to the w all have a seventh-power profile. That is. 


for 0 < < 1, 


u ~ "w 
"e - "w 


r 1 / 7 


S T 


for 1 < £i 


(40) 


*E> 


u - 




u e " "w 


= 1 


J 


It is also necessary to express T/T e (= p g /p) as a function of £g,. 
The; form assumed herein is 


for O^g^l, 




± !>L w 1 / 7 ,t z /A 

p - T T e + b ^T " c£ I ) 


for l<fc 


T J 


where 


Pe T 

« 77T“ « 1 

P T e 


T 

T 

.--- 1 


(41) 
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c 



*) 


T, 



For a Prandtl number of 1, the dependence of T/T e on u/u e in equations 

( 41 ) is the same as that indicated in equation (l2). When the Prandtl 
number is not too far from 1, equations (4l) should give at least a 
reasonable estimate for T/T g . 
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By using equations (40 ) and (4l), the boundary-layer momentum thick- 
and displacem e nt thickness can he expressed as 



(42a) 


(42b) 


where I g , Iy, and I Q are functions of h and c defined by the 
integral 


I N 


r 1 n, 

/ z dz 

Jn 1 + bz - 


cz 


(43) 


with N« 6, 7, or 8. A method for evaluating ^ is indicated in 
appendix D. It should be noted that the integrals ly and (ly - I Q ) 
are tabulated in reference 11 for the range -1 < b < 10 and 

0 < ^ 1 j - < 1. The results were computed using Simpson's rule, re- 

sults of reference 11 can be vised to tabulate values of I 8 ftnfl I 

The integral I g is found from the relation I g » y - bI 7 + cI 8^ wh ich 

follows from equation (43) . The reciprocals of I g , I ? , and I g , based 
on the data of reference 11, are presented in table IV. Reciprocals axe 
used to permit linear interpolation. When b is near -1 or 
c/(b + l) is near 1, the linear interpolation becomes inaccurate and use 
of appendix D is recommended. 


Since 0/5 is independent of x (assuming the wall temperature to 
be constant), equation (39) can be written 


“w 


P u 2 
H e e 


0 dS 
5 


(44) 


In order to integrate equation (44), it is necessary to have a relation 
between t ¥ and 6. For incompressible turbulent flow past a semi- 
infinite flat plate, the Blasius relation (ref. 12) is 


*w 


P u 
e e 



( 45 ) 
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This fo rm ula can he extended to approximate the compressible flow over 
a semi-infinite plate by evaluating the fluid properties at a suitable 
mean temperature T m (refs. 10, 11, and 13) giving 

1/4 

~~~~S “ 0.0225 (46a) 

p u* \u By 

me v e / 


or 


where 



(46b) 


From reference 13, a reasonable estimate for T is 
' m 


T m - 0.5(T W + T e ) + 0.22(T r - T 0 ) 


(47) 


When the wall is moving, the logical extension of equation (46b) is to 
use velocities relative to the wall. Thus, 


or 



(48a) 



(48b) 


Equations (48) are expressed, as indicated to ensure that has the 

same sign as u g - u^., which is to be expected from physical considera- 
tions. Substitution of equation (48b) into equation (44) permits the 
latter to be integrated, for T w constant, giving 
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Other boundary-layer parameters of interest are 


'W 


0.0460 


p a u 
^e e 


(50a) 


e _ dS 
u e “ dx 


0.0460 ^<p 


- V^ 4/s 


e/S 


- V’ 


■*-' fe) 


1/5 


(50b) 


The heat transfer at the wall can be estimated using a Reynolds 
analogy. For compressible flow over a semi- infinite flat plate, the 
Reynolds analogy may be expressed as (ref. 13) 


«*■ 


S.mC^V - T r> „ -«/* 


u 


e 


m 


w 


(51) 


When the wall moves, an estimate of the heat transfer might be obtained 
from 




c (T - T )o 
p,m N w r / m 

u e " ”w 


■2/3 


w 


(52) 


It is pointed out in the CORRELATION OF DATA, section that the exponent 
of cr^ should probably decrease in magnitude with increasing u^/u e . 

The recovery temperature T is yet to be determined. By analogy 

with the turbulent flow over a semi-infinite flat plate, T r may be 
estimated as (ref. 13) 


T r (** A 


2 Ugr(0) 

zFc"~ 

e 


(53a) 


where 


r(o) s 




(53b) 


Equation (53a) is solved by iteration since c^ m and r(o) are func- 
tions of T r . For the laminar -boundary-layer case, r(o) increased with 
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an increase in u^Ug. This suggests that equation (53b) might under- 
estimate r(o) at the higher values of u^/ug. However, equation (53b) 

is probably -well within the accuracy of the other assumptions required 
to solve the turbulent -boundary-layer equations. 


j surface temperature . - For convenience, it was as sinned during 
the course of the turbulent -boundary-layer solution that the wall sur- 
face temperature was constant with x. The validity of this assumption 
can be investigated by solving for the wall temperature as was done for 
the laminar boundary layer. 


Equating equations (52) and (30), plus equation (48b), gives 


*w- T b 

*b 


BtCnyog - i] 

1 + B 


= B 



+ 0(B 2 ) 


(54) 


where 


B = 0.0815 



V U € 


.3/5 




3/10 


TTnUTre the laminar-boundary-layer case, B varies with x. Therefore, 
T/O^ varies with x, which contradicts the original assumptions under 

which equations (54) were derived. However, equations (54) can be used to 
establish the conditions under which T^ is nearly equal to T^. The 

procedure to follow is to evaluate the right side of equations (54) 
assuming T w = 3^. (Eqs. (33) are applicable.) If the resulting value 


for (0^ - is small (e.g., less than 0.1 for all values of x 

under c onsider ation), then the assumption of a constant wall temperature 
is reasonably accurate and the solution indicated by equations (54) is 


valid. 


The conditions under which the wall temperature behind a strong 

shock wave (M 2 » 1, u^/u -*• 6) remains essentially at its original 
value will now be determined. Assume c^ is constant, a - 1, y *» 1.4, 
p. ia proportional to T, T^. - T^, and 1^ =* l/£(l + b - c) (N + l)J. 

The latter is the leading term in the expansion for Igj when equations 

(D3) and (D4) of appendix D are used. The following relations can then 
be obtained: 
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9/8 

1 - V u e 


7 + 2(V u e) 




72 


7§ for V u e 


Tm 5 - 56 K/ u e ) - 0.28 [1 + (u,j/ u e )^] 

r e ecys,) “ 1 



for (u^/ug) a 6 


( 

*“Vv 


1/2 


f 


for (u/u e ) a 6 


(55) 


Evaluating B for il/u » 6 and substituting into equations (54) yield, 
for^»l. 


^w “ T b g 0.08 


*b 


Assume the fluid to be air and the wall to be cast iron, both initially 
at I, a 530° R. it is consistent with the assumptions of p oc I, a a l, 

. 15: V 

and so forth, to consider ^ AJ— — ■£ — *» 0.0004, 

V w 0" 

using the property values of reference 7 for 3^ « 3 ^ a 530° R. Equa- 
tion (56) becomes 

(57) 


Thus, when is large, the wall surface temperature will rema i n within 

10 percent of its original value for (xu e / v e ) 3 / 10 l^ <0(l0 4: ). 

( T r he preceding solution requires that the wall thickness be at least 
greater t han the wall thermal -boundary-layer thickness (eq. (36)). 
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■ CORRELATION OP DATA 

In order to correlate the data pertaining to the boundary layer 
behind a wave, it is necessary to define a Reynolds number which char- 
acterizes the flow. 

Consider the t ,x-diagram for the flow behind a wave (sketch (a)). 



An expansion wave is indicated in sketch (a) for convenience. When a 
particle of the external stream moved a dis tance K » (u e - u^Jt 

relative to the wall, it is a distance x = £ + u^t, or x » £/(l - u^/ug), 
behind the wave. Thus, K ~ x(l - u^/ug) is a physical a-lm^rmion (corre- 
sponding to a given x) which characterizes the boundary-layer develop- 
ment. Similarly, (u^ - u g ) the velocity of the external flow relative 

to the wall is the characteristic velocity. The Reynolds number may 
then be defined as Re a ?C u e ~ or 

Re = (u e x/v w )(l - u w /u e ) 2 (58) 

where v is arbitrarily referenced to the wall temperature. Equation 
(58) reduces to Re = (u e x/v w ) for u^ » 0. Equation (58) differs by 

the factor Ug/u^ from the Reynolds number defined by equations (32) of 
reference 1. For u e = 1, the Reynolds number of reference 1 agrees 
with equation (58), but for other values of u^u e , the two Reynolds 

numbers differ. In particular, the Reynolds number of reference 1 be- 
comes infinite as u^ -»■ 0. Hence, equation (58) appears to be a more 
useful form. 

I 

T he correlation of the laminar-boundary-layer skin-friction and 
heat-transfer data will now be indicated. The local s kin -friction coef- 
ficient may be defined as 


C45 ~ 




f " 1 , ,2 
2 p vK “ V 


( 59 ) 
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Therefore, using equations (l6a), (58), and (59), 


,-y^Re 


-Jz f'Co) 
IV u e - 1 | 


(60) 


where c^-^Re is positive if the wall shear is in the +x-direction. 
Equation (60) is positive for 0 < '%/ VL e < 1 and. negative for 
u^/ug > 1. The magnitude of c^-^/Re increases monotonically from 0.664 
at u^/Ug b 0 to 1.128 at u^/ug « 1 and 2.291 at u^/ug « 6. Using 
equations (16) and (27), the Reynolds analogy can he written as, 

for 0 < u^/ug < 1, 


^w 


- u (a ) 0 - 6s -°- 15 (v' u e) 
e v w' 


T w T w “ T r °p,w 


B 1 


(61a) 


for 1 < u^/ug < 6, 


-v - % 


*c T - T c 
w w r p,w 


» 1 


(61b) 


The exponent of cr^ is seen to decrease from a value of 0.65 at 
u^/Ug = 0 to a value of 0.39 at u^/ug « 6. 

By using equation (50a), together with the definitions of Re and 
c f as given by equations (58) and (59), the turbulent skin friction 

can be expressed as 


CjRe 


l/ 5 


0.0920 


Aw0 1 Vm/PnA 4 

‘ 1 - 

"e 

l T e 6-! . Sr V P J 


i.i 

-V u J -1 


u e 


(62) 


The Reynolds analogy proposed for the turbulent case is given by equa- 
tion (52), which may be rewritten as 


Sf K - %) ° m 2 / 3 

T w ^ " T r) e p,m 


« 1 


(63) 
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Comparison of equations (6l) and (63) suggests (except for the use of a 
mean reference temperature in eq. (63)) that the exponent of c m should 

decrease with increasing u^/u g . However, the use of the exponent 2/3 

is probably within the accuracy of the other assumptions involved in the 
turbulent -boundary-layer analysis. 


In reference 4, the boundary layer behind a wave was estimated by 
using the concept of "an equivalent semi-infinite flat plate." That is, 
reference 4 assumed "... the properties of a lamina of fluid in the shock 
tube which has been in motion with a velocity (u g - Uy.) for a time t 

are equivalent to those of a lamina of fluid which has progressed rear- 
ward for a period of time t from the leading edge of a semi-infinite 
flat plate in a steady flow with free -stream velocity (u 0 - u^) . . . . " 

For the laminar case, this is equivalent to replacing equation (60) by 


Cp-^/Se = 0.664 


I - %/ u e 

I I - %/"e| 


(64) 


whi* is accurate for « 1, but underestimates the shear by a 

factor of 1.7 at U v /' li e ~ 1 813(1 "by a factor of 3.4 at u^/iig m 6. For 

the turbulent case, this is equivalent to using equation (62), but with 

\/<i ■ 53 evaluated for u^/ug » 0. If the variation of fluid state 

properties is neglected, equation (62) can be expressed as 

\l/5 1 




1/5 


0.0577 


MS) 


- (W 

|1- K/-e)| 


(65) 


which has the values 0.0577 , 0.0607, and 0.0705 at x \ l / u e 83 0, 1, and 6, 

respectively. The fact that the values of c^Re^/^ for u^/ug » 0 and 

1 differ relatively little suggests that the equivalent flat-plate ap- 
proach may give reasonable estimates for the turbulent boundary layer be- 
hind a weak wave. Another interpretation is that (by comparison with the 
laminar solution) the present method underestimates the wall shear for 
values of Uy/ug other than those satisfying u^/u fe « 1* Note that 

the seventh-power velocity profile and the Blasius relation between 
shear and boundary-layer thickness (eq. (45)) apply to fully developed 
turbulent pipe flow as well as to the turbulent flow over a semi-infinite 
flat plate. It is hoped that these expressions are sufficiently univer- 
sal so that they can be used herein without modification, except to use 
velocities relative to the wall. The validity of this approach must 
ult imat ely be verified by experiment. 
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CONCLUDING REMARKS 

The bou ndar y layer behind, a shock or thin expansion wave moving 
into a stationary fluid, has been studied.. Both 1 ami mr turbulent 
boundary layers were considered. The wall temperature behind the waves 
was also Investigated. 

For strong shock waves, the boundary layer behind the shock will 
generally have very large temperature gradients normal to the wall. The 
assumptions upon which the laminar-boundary-layer solutions were obtained 
(eqs. (8)) may then have to be reevaluated. It may be possible to obtain 
more accurate estimates of the i aminR-r boundary layer by uping some wie>nn 
reference temperatures as discussed in appendix C. The choice of the 
appropriate mean temperature might be obtained from either experimental 
data or by comparison with more accurate integrations of the laminar- 
boundary-layer equations. 

The turbulent-boundary-layer solutions presented herein represent 
an extension of empirical , semi -Infinite flat-plate, boundary-layer data 
to the case where the wall is moving. The basic assumptions (especially 
those incorporated in eqs. (40 ), (4l), and (48)) mus t ultimately be ver- 
ified by experiment, particularly for the case of strong shocks . 

It should also be noted that the fluid behind a very strong phpr»k 
may be dissociated. Such dissociation introduces complexities of the 
boundary-layer problem that are beyond the scope of this report . Exper- 
imental studies of the boundary layer behind very strong shocks would 
yield fundamental information relating to boundary-layer development in 
a dissociated fluid. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, February 16, 1956 
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APPENDIX A 
SYMBOLS 

function defined by eq. (32b) 
local speed of sound 
function defined by eqs. (54) 
function defined by eqs . (4l) 
function defined by eqs . (41 ) 

/l 2 

local skin-friction coefficient, t /— P w (u w - u g ) 
specific beat at constant pressure 
complimentary error function 
function of tj defined by eq. (6) 
integral defined by eq. (43) where 11=6,7, or 8 


thermal conductivity 


Mach number in x,y-coordinate system, M g = (u e /a e ), etc. 
pressure 

rate of heat transfer in +y-directian 
gas constant 


Reynolds number 


’ % v u . 


s ) 


function defined by eqs. (15) 
function defined by eqs . (15) 
absolute static temperature, °R 

absolute mean static temperature (e.g., eq. (47)), °R 


absolute wall surface static temperature for zero heat transfer, 
°R 


t 


time 
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u,v 

u,v 

u d 

x >y 

x,y 

a 

r 

Ap 

6 

Ap 

6 

8* 

s 

£n 


3 T 


^8 


velocities parallel to x,y-axes 
velocity of wall In x, y-coordinat e system 

velocities parallel to x,y-axes 

velocity of wave in x,y-coordinate system 

coordinates stationary with respect to wave (fig. l(b) ) 

coordinates stationary with respect to wall (fig. 1(a)) 

diffusivity, k/pc^ 

ratio of specific heats 

wall thermal -houndary-layer t hickn ess corresponding to 


- T w 


= 0.99 


fluid velocity-houndary-layer thickness corresponding to 
u - iv 


u e ‘ % 


= 0.99 


fluid thermal-houndary-layer thickness corresponding to 


T - T. 


w 


T - 
e w 


= 0.99 


fluid velocity-houndary-layer thickness from Pohlhausen viewpoint 
fluid houndary-layer displacement thickness, / { 1 \ dy 




similarity parameter for integral solution of laminar case, tj/t) 
similarity parameter for integral solution of turbulent case, 

y/8 

similarity parameter, eq. (5) 
value of T) corresponding to y = 8 

fluid boundary -layer momentum thickness, I pu ( 1 - — J dy 

\J 0 p e u e V u e/ 


<60 
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p coefficient of viscosity 

v kinematic viscosity 

p mass density 

a Prandtl number, pc^/k 

local shear stress exerted by fluid on vail 
(p function defined by eq. (46b) 

ijr stream function (eqs. (4)) 

Subscripts: 

b undisturbed flow ahead of wave 
e flow external to fluid boundary layer 

m quantity evaluated at 

w quantity evaluated just above wall surface, y = 0 + 

0“ quantity evaluated just below wall surface, y = 0~ 

Superscripts: 

’ denotes different iat ion with respect to tj 

Special notation: 

0( ) order of magnitude 
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APPENDIX B 


WAVE RELATIONS 

For convenience, useful equations relating conditions across a wave 
are noted herein. Shock waves and expansion waves are cons ider ed 
separately. 

S Shock wave (u^/ug > l) . - Consider the flow in the x,y-coordinate 

10 system, (.i-e., coordinate system moving with the shock) . Let subscript 

h designate the undisturbed flow upstream of the shock wave »riH 


* 


7, 7 


£ 

II 


u e 

V *b> 


T e> a e* Pe 





v / / / / / yv//\ 1 


x, u 


subscript e the flow downstream of the shock wave and external to the 
boundary layer. Note that u^ = u^. so that = u^/^ = u y/%- Then, 
from ordinary shock-wave theory. 


(r + i)m| 

u e Cr - i)m^ + 2 




> 




+ 5 


for r = 1.4 




<- 


2u w/ u e 


Cr + 1) - (r - lJCu^/ug) 

_ S V u e 


6 - (V u e> 


for r = 1.4 


J 


(Bl) 


(B2) 
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M 2 = 




e Cr + lJCuyug) - (r - i) 


for x 


T, 


6 K/ U e ) - 1 

(r + lJCu^/ug) - (r - i) 




J 




( u w / u e )[fr + l) - (r - lKu^/ug)] 

6 ( u w/ u e) " 1 


> 


CV u e >[ 6 - 


for x — 1.4 


(B3) 


(B4) 


Expansion wave (u^/u^) < l) . - The expansion wave is assumed to have 

negligible thickness. In the x,y-coordinate system (i.e., fixed with 
respect to wall), the actual and assumed form of the expansion wave is 
as follows : 


Expansion 

wave 


u d B 

V Pb 


u 


T e J a e> P e 


'/// s/y/ 7 

Actual expansion wave 


-“b 

u e 

V V *b 

T e> a e j P e 

//////// 

V' ////// 



Assumed expansion wave 


The discontinuity moves with velocity = -a^ into the undisturbed 

fluid. Conditions on either side of the assumed expansion wave are the 
same as shown for the actual wave. In the x,y- coordinate system the 
following figure applies: 



3959 



CT^ "back 


• MCA TN 3712 


35 


The state properties across the wave are defined by the isentropic 

r-i r-i 

relations a e /a b = (Tg/0^) 1 / 2 = 2 = ^ • Other equa- 

tions of Interest are 


Pe 

Pb 


1 T - 1 1 I K/ U e>' 

2 V u e 


- l] ? 
L 5 V u e J 


2r/(r-i) 




for x ~ 1.4 






^ (r - i) 

(r + i) - 




6 " 5 (P e /Pb> 


iT? 


for r = 1.4 




J 


M = 
e 


(r + lJCuy/ug) - (r - i) 

5 


> 


6 K/ U e) ■ 1 


for x = 1.4 


J 




5»_ 1 r - 1 1 - K/ U e> 

®b ■ ’ 


2 “7% 

2 


6 C u w/ u e ) - 1 

L ®V u e I 


for x - 1.4 


(B5) 


(B6) 


(B7) 


(B8) 
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APPENDIX C 


LAMINAR -BOUNDARY -LA2ER SOLUTION REFERENCED 
TO A MEAN TEMPERATURE 

In the body of this report, the laminar-boundary -layer solution 
was obtained, by assuming that p. and. k are proportional to T and. 
that a and. c are constant, all referenced, to wall conditions 

p 

(eqs. (8)). The question naturally arises as to whether the -use of a 
mean reference temperature ^ would lead to more accurate correla- 
tions with exant (in regard to fluid properties) integrations of the 
laminar -boundary-layer equations. 

If a mean reference temperature is vised in equations (l6), they 
become 


: v - "(°> || 


^PmFm 

2x 




r dr] + 




s d^ 


where 


T (m^ u|r(0) 

r = 1 + (« ^ 2tV“ 

e \ e / e p,m 


reference 13, it is found that the choice 


T m = 0.5(T w + T e ) + 0.22(T r - T g ) 


(Cl) 


(C2) 


( 05 ) 


(04) 


and the solutions for r and s are considered functions of o m . In 


(05) 
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correlated equations (Cl) and (C2) with accurate integrations of the 
"boundary -layer equations. Plow past a semi- infinite flat plate is 
considered therein. Equation (C5) should also provide improved estimates 
for and q^ for u^/ug ^,0. If v e /u e is to "be correlated with 

accurate integrations of the "boundary-layer equations, another estimate 
for mi ght "be called for. Note that equation (C5) defines the 

temperature which was used in the turbulent-houndary-layer analysis. 


o 
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APPENDIX D 


EVALUATION OP 


% 


The turbulent -houndary-layer analysis requires evaluation of the 
integral 


, r 1 As 

vJq 1 + hz - cz 4 


•where N = 6, 7, or 8, and 
Q 

t= ^" 1 

The denominator of equation (43) may he factored, so that 

o 

(-c)(x - X)(x - p) = 1 + hx - ex where 


(43) 


X 

P 


- i (* + 4c ) 

=s( b -'f 2+4 =) 


from physical arguments h > -1, c ^ 0. and 0 ^ y- c -■■■ < 1. It can then 

1 + b 

he shown that 0 < 0 and X > 1. Equation (43) can now he written as 






+ 4c 


^P " 


(Dl) 


where 


T _ , z^dz 
T^P ~ z - P 
'0 


p 1 N 

^ m Jo - 
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o 


3!he integrals and can be integrated as f oIIowb : 


■-'[-M-ekb] 

m=l -* 


dz 


(D2a) 


Similarly, for 



(D2b) 


The logarithm term is very nearly equal and opposite to the summation 
term, in equations (D2a) and (D2b), for values of p not close to zero 
and values of X not close to 1, respectively. Hence, each of these 
terms must be evaluated very accurately in order to get accurate esti- 
mates of and Ijj For values of p not close to zero nnrt val- 

ues of X not close to 1, alternate expressions for I w R and I w t, 
are desirable. Since p < 0, * 


- 


1_ P 1 z H dz 

"I RH 


r r 


- P 


m=0 '-'0 


..... 1 V N.'ml f 1 

i-pL (n + m + i; i vi - p )’ 

m=0 

= i r + i • 

W + IKI - pj|_ (n + 2 hi - p; 


21 


(JST + 3)(N + 2)(l - p)' 


] 


(D3) 
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Since X > 1, 




dz I 


w 


w 


co 

+ DU - X) |l - + 1) Ys TF+'m)(N + m + 1) (£) J 
i r (n + 1) 

+ l)(l - X) (n + l)(N H 


. (N + 1) 

2 ^ X (N + 2)(N + 3)X 2 

(D4) 

Equations (D3) and (D4=) converge rapidly for p « 0 and X » 1. 


3 


REFERENCES 

1. Mir els, Harold: Laminar Boundary Layer Behind Shock Advancing into 

Stationary Fluid. HACA OUST 3401, 1955. 

2. Donaldson, Coleman duP., and Sullivan, Roger D.: The Effect of Wall 

Friction on the Strength of Shock Eaves in Tubes and Hydraulic 
Jumps in Channels . BACA TUT 1942, 1949. 

3. Holly er, Robert B., Jr.: A Study of Attenuation in the Shock Tube. 

Eng. Res. Inst., Univ. of Michigan, July 1, 1953. (U.S. Navy Dept., 

Office Naval Res. Contract No. N6-0NR-232-T0 17, Proj . M720-4. ) 

4. Trimpi, Robert L., and Cohen, Nathaniel B.: A Theory for Predicting 

the Flow of Real Gases in Shock Tubes with Experimental Verification. 
HACA TN 3375, 1955. 

5. Chapman, Dean R., and Rubes in, Morris W. : Temperature and Velocity 

Profiles in the Compressible Laminar Boundary Layer with Arbitrary 
Distribution of Surface Temperature . Jour. Aero. Sci., vol. 16, 
no. 9, Sept. 1949, pp. 547-565. 

6. Low, George M. : The Compressible Laminar Boundary Layer with Heat 

Transfer and Small Pressure Gradient. BACA TN 3028, 1953 

7. Hartunian, Richard A.: Some Viscous Effects in Shock-Tube Flow. 

M- S. Thesis, Cornell Univ.,- June 1954. 


3959 



CT-6 


MCA TN 3712 


41 


8. Rott, N., and Hartunian, R.: On the Heat Transfer to the Walls of a 

Shock Tube. Graduate School of Aero. Eng., Cornell Univ., Nov. 
1955. (Contract AF 33 ( 038) -21406. ) 

9. Eckert, E. R. G.: Introduction to the Transfer of Heat and Mass. 

McGraw-Hill Book Co., Inc., 1950, pp. 266-275. 

10. Tucker, Maurice: Approximate Calculation of Turbulent Boundary- 

Layer Development in Compressible Flow. MCA TN 2337, 1951. 

11- Bartz, D. R.: An Approximate Solution of Compressible Turbulent 

Boundary -Layer Development and Convective Heat Transfer in 
Convergent -Divergent Nozzles. Trans. A.S.M.E., vol. 77, no. 8, 
Nov. 1955, pp. 1235-1245. 

12. Sc hli chting, H. : Lecture Series "Boundary Layer Theory." Pt. II - 

Turbulent Flows. MCA TM 1218, 1949. 

13. Eckert, Ernst R. G.: Survey on Heat Transfer at High Speeds. Tech. 

Rep. 54-70, Aero. Res. Lab., Wright Air Dev. Center, Wright- 
Pattersan Air Force Base, Apr. 1954. 



42 


mck ts 3712 


TABLE I. 


LAMINAR BOUNDARY LAYER BEHIND WAVE 
(a) Solution for Uy/Ug => 0 


(EQS . 9 AND 15) 
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TABLE I. - Continued. LAMINAR BOUNDARY LAYER BEHIND WAVE (EQS . 9 AND 15) 

(c) Solution for u^/Ug “ 0.50 
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TABLE I. - Continued. LAMINAR BOUNDARY LAYER BEHIND WAVE (EQS . 9 AND 15) 


(d) Solution for u^/Ug =0.75 
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.00 0 0 

14 96 3 

.0 005 

.0016 

14)038 

4)003 

.0 010 

<5 

4.8 05 5 

1.0000 

.00 0 0 

14 9 6 3 

.000 4 

.0012 

14)038 

4)003 

.0 00 8 


4.3 8 55 

10000 

4)000 

14. 96 4 

.0 00 3 

.0008 

14)038 

.0 0 02 

.0 00 6 

n 

mil 

£8888 

$888 

£1113 

:888i 

$881 

£8838 

$881 

$881 

4J9 

4.6 85 S 

10000 

4)000 

1.1964 

.0001 

.0003 

14)038 

4)001 

.000 2 

8.0 

4.7 85 5 

10000 

.00 0 0 

11964 

.0001 

.0 00 3 

14)038 

.OOOO 

.0 00 8 

S.1 

4.8 05 5 

10000 

.0000 

1196 4 

.000 0 

.00 0 3 

UO 0 3 8 

.0 00 0 

.0 00 1 

If 

8IH 

£8888 

$888 

am 

.0000 

.0600 

.00 01 
.0 0 01 

£8831 

.0000 

.OOOO 

.0 00 1 
.0 00 1 

5.4 

5.1055 

10000 

4)000 

11964 

.0000 

.0001 

UO 03 8 

.0 00 0 

.0 00 0 

SJ5 

5.2 05 5 

1.0000 

.0000 

1196 4 

.0000 

.000b 

14) 03 8 

.0 00 0 

•000 0 
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TABLE I. - Continued. LAMINAR BOUNDARY LAYER BEHIND WAVE (EQS. 9 AND 15) 


(e) Solution for %/ug = 2 (ref. 1) 






| Prandtl number v M , o.72 

1 

f 

f* 

f" 

HB 

r 

-r ' 

psa 

8 

n 

0.0 

0*0 0 0 0 

2.000 0 

-13 19 1 

0.0 00 0 

0.8 9 97 

0.0 0 0 0 

0.0000 

3.0000 

0.8512 

•1 

.19 4 9 

1.8 9 84 

— 3 jO 0 9 1 

.089 7 

5922 

.14 7 9 

.0958 

.9151 

.8 4 52 

.2 

J790 

3-7 9 8 8 

“5 8 04 

.178 0 

.8 7 04 

.2 861 

.1831 

.8 313 

.8 27 9 

3 

5548 

3-7 0 29 

-5 35 6 

.2634 

.8 356 

.4 0 66 

.2621 

.7 4 98 

.8004 

A 

.7 20 5 

1.6 1 2 1 

-5777 

.344 7 

.7898 

5043 

.3331 

.6715 

.7645 

JS 

5774 

15 27 6 

-5103 

.4210 

.7356 

5758 

.3965 

.5973 

.7217 

.6 

1.0 26 3 

1.4 503 

-.7 367 

.491 6 

.6 755 

5 309 

.4527 

.5374 

.6739 

.7 

U. 677 

15804 

-j6 6 01 

.556 0 

.6122 

.6 413 

.5022 

.4 6 25 

.6227 

JB 

13036 

1.3183 

-58 34 

.614 1 

•5480 

5398 

.5454 

.4029 

.5 69 7 

S 

3-4 316 

1.2 63 7 

-5088 

.665 7 

•4849 

.6206 

.5839 

.3 4 86 

.5162 

3LO 

15556 

i.ai64 

-.4 3 8 2 

.7111 

.4 24 4 

.5878 

.6153 

.2996 

.4 63 6 

u. 

16751 

1JL 7 5 9 

-5728 

.750 7 

5676 

.5 4 55 

.6430 

.3 5 58 

.4127 

Z.3 

3-7910 

3-1416 

-5135 

.78 4 8 

.315 5 

.4 97 4 

.6666 

.2170 

.364 3 

13 

13036 

1.1 12 9 

-5606 

.8139 

.2683 

.4 4 65 

.6866 

.1829 

.318 9 

1.4 

3.0137 

3-0 8 93 

-5142 

.8 38 6 

J3262 

.3953 

.7034 

.1531 

.876 9 

13 

2.1 21 6 

1.0 6 99 

-1742 

.8 59 4 

4. 8 93 

5456 

.7174 

.127 4 

.23 8 6 

13 

2527 8 

15 54 2 

-1403 

.876 6 

.15 70 

.2988 

.7290 

.1053 

.30 4 0 

1.7 

33336 

1.0 417 

-1116 

.89 09 

.1293 

J3 556 

.7385 

.0 8 64 

.1731 

13 

2.4 36 2 

1.0 317 

-JD 879 

.9 02 6 

.1057 

5166 

.7463 

.0 705 

.145 8 

13 

253 9 0 

1.0 2 3 9 

-5686 

.912 2 

.0 6 58 

.1819 

.7527 

.0572 

.1219 

33 

33411 

1.0 1 7 9 

-.0 5 2 9 

.9 19 9 

.0692 

.1514 

.7578 

.0 4 60 

.1013 

34. 

2.7 4 2 6 

1.0 1 3 3 

-0404 

.9 261 

.0 55 4 

.1250 

.7630 

.0 3 68 

.08 3 3 

33 

25 438 

1.0 0 9 7 

-.0 3 0 6 

.9 310 

.0 4 41 

4.024 

.7653 

.0 29 3 

.0 68 2 

33 

25 4 46 

1.0 0 7 1 

-0229 

.935 0 

.0 348 

.0 832 

.7679 

.0 2 31 

.055 3 

3.4 

3.0 4 52 

1.0 0 51 

-.017 0 

.9 38 1 

.0 2 73 

.0 671 

.7699 

.0181 

.04 4 6 

33 

34 4 5 6 

13036 

-.0 12 4 

.9 4 05 

3 213 

.0538 

.7715 

.0141 

.0357 

33 

33 4 59 

1.0 026 

T 1 1 

.94 2 4 

.0165 

.0 4 27 

.7738 

.0109 

.0284 

3.7 

33 4 62 

1.0 0 1 8 


.94 3 8 

.0127 

.0337 

.7737 

.0 0 84 

.0 23 4 

33 

33 463 

13012 

-50 4 6 

.9 44 9 

3097 

.0 264 

.7745 

.0064 

.017 5 

33 

33464 

150 0 8 

-3 0 33 

.9 45 8 

.0073 

.0 20 5 

.7750 

.0 0 49 

.0136 

33 

33465 

13006 

“5 0 23 

.9464 

JO 055 

.0158 

.7755 

.0037 

.0105 

34 

3.7 4 65 

15004 

-3 016 

.946 9 

3041 

.0121 

.7758 

.0 0 28 

.0081 

33 

33466 

15 0 0 3 

-3 011 

.9473 

3031 

.0 09 3 

.7760 

.0020 

.00 61 

33 

39 466 

15 0 0 2 

-3 0 07 

.9475 

3023 

.007 0 

.7763 

.0015 

.00 4 6 

33 

43466 

14)001 

-3 0 05 

.9477 

3017 

.0 0 52 

.7763 

.0 011 

.0035 

3*5 

44466 

13001 

-3 0 03 

.9479 

3012 

.0 0 39 

.7764 

.0008 

.0026 

33 

45 4 66 

13000 

-jO 0 0 2 

.9 48 0 

.0 0 09 

.0 0 29 

.7765 

.0006 

.0019 

3.7 

43466 

15 00 0 

-3 0 01 

.9 48 0 

.00 06 

.00 21 

.7765 

.0 0 04 

.0014 

33 

43466 

15000 

-,0 0 01 

.9 4 81 

.0 0 05 

.0 015 

.7766 

.0 0 03 

.0010 

33 

45 4 66 

1.0 0 0 0 

-3 0 01 

.9481 

.0003 

.0011 

.7766 

.0 0 02 

.000 7 

4.0 

45466 

15 00 0 

4) 0 00 

.9 48 2 


.00 08 

.7766 

.0001 

.0005 

44 




.9 48 2 

BTTTrrH 

.00 0 6 

.7766 

.0 0 01 

.0004 

4.2 




.94 8 2 

.0 0 01 

.00 04 

.7766 

.0001 

.000 3 

43 




.9 48 2 

.0 0 01 

.0 0 03 

.7766 

.0001 

.0002 

4.4 




•9 48 2 

.0 0 01 

.0 0 02 

.7766 

.0000 

.00 01 

45 




.9 48 2 

.0000 

•0001 

.7767 

.0000 

.0001 

43 




.94 8 2 

4)000 

.0001 

.7767 

.0000 

.0001 

4.7 




.9 48 2 

.0000 

*0 0 01 

.7767 

.0000 

.000 0 

45 

45 




.948 2 

4)000 

.00 0 0 

.7767 

.0000 

.0000 
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TABLE I. 


Continued. LAMINAR BOUNDARY LAYER BEHIND SHOCK WAVE (EQS. 9 AND 15) 


(f) Solution for u„/u e = 4 (ref. l) 






Prandtl number o v , 0.72 

n 

f 

f ■ 

f 

J7 rdn 

r 

-r* 

X' 8<w 

s 

-»» 

0.0 

0.0 00 0 

4D000 

-40 6 2 3 

0.0 00 0 

0.9129 

0.0000 

0*00 0 0 

10000 

14.156 

.1 

0798 

30964 

-35 8 45 

.0 9 09 

.8996 

.2582 

.09 45 

.8 892 

14.005 

.2 

.719 8 

32077 

-3.7 7 02 

.1791 

.8623 

.4 847 

.1779 

.7811 

ID 5 7 6 

.3 

to 3 3 8 

30 457 

-34 5 46 

.262 6 

.8047 

.6579 

.2509 

.6 78 4 

3 9 31 

.4 

U3901 

36188 

-30 76 6 

.339 6 

.7328 

.7690 

-3139 

•5 6 30 

3 136 

* 5 

10373 

32 313 

-20 717 

..408 9 

.6529 

.8202 

.3677 

.4960 

.8253 

.6 

1.7 377 

1084 4 

-22 6 88 

.4701 

.5705 

.8209 

.4134 

.4180 

.73 3 7 

•1 

10355 

1.7 7 6 8 

-1.8 8 8 7 

.5231 

.4900 

.7838 

.4517 

.3 4 92 

.64 39 

.8 

2.0 9 4 3 

100 5 5 

-15 4 4 6 

.5683 

.4 145 

.7318 

.4 8 35 

.3 893 

.5 56 3 

.9 

2.2 4 7 6 

10664 

-10 430 

.606 2 

.3461 

.6458 

.5098 

.2 378 

.4 75 7 

duo 

30885 

U 5 5 4 

-5858 

.637 8 

.3855 

.5646 

.5313 

.1 939 

.4036 

JL1 

80195 

10679 

-.7 713 

.663 6 

.3331 

.4 844 

.54 8 8 

•15 7 0 

.3374 

1.8 

2.6 4 27 

1J. 9 98 

-.5 9 5 7 

.684 6 

.1885 

.4089 

.56 3 9 

126a 

.3801 

1.3 

3.7 6 0 0 

13 47 6 

-A 547 

.7016 

.1511 

.3406 

.57 4 3 

1007 

.2306 

3.4 

20737 

1J. 07 9 

-3431 

.7151 

.1201 

.2 804 

.5833 

.0 79 8 

.1883 

1*5 

3.9819 

ID 7 81 

-2 56 0 

.735 8 

.0948 

.2283 

.5903 

.0 6 28 

3.535 

1.6 

3D 8 8 5 

ID 5 6 0 

-.1 8 9 0 

.734 3 

.0 742 

.1842 

.5959 

.0 4 91 

3226 

1.7 

33 9 3 3 

1.0 3 9 8 

-4.38 0 

.740 7 

.0 5 77 

.1473 

.6003 

.0 38 2 

.0 97 8 

1.8 

3D 9 67 

10 38 0 

-2 99 8 

.745 8 

.0445 

.1168 

.6036 

.0 2 94 

jO 77 4 

3L9 

30990 

IX) 19 5 

-0714 

.749 7 

.0341 

.0919 

.6 0 6 2- 

JO 2 26 

.0608 

24) 

36006 

10135 

-.0 5 06 

.753 7 

.0 260 

.0717 

.60 6 2 

.0172 

4)474 

23 

3.6018 

1.0 0 9 3 

-.0 355 

.755 0 

.0196 

.0556 

.609 7 

.0130 

D 3 6 7 

2J3 

3.7 0 2 5 

10 06 3 

-.0 24 6 

.756 7 

.0147 

.0437 

.610 8 

.0 0 97 

4)282 

80 

30030 

1.0 0 4 3 

-0 169 

.758 0 

.0110 

.0326 

.6116 

.0073 

D216 

2*4 

3.9 0 3 4 

10038 

-0115 

.758 9 

.0083 

.0347 

.6133 

D 0 54 

4)163 

2^ 

44) 0 3 6 

10016 

-0077 

.7596 

.0060 

.0186 

.612 7 

.0040 

4)123 

3.6 

4.103 8 

ID 013 

-6053 

.7602 

.004 4 

.0139 

.6131 

JO 0 29 

D 0 9 3 

2.7 

42 0 39 

1J0 0 08 

-0034 

.7605 

.0033 

.0103 

.6133 

D 0 21 

JO 0 6 8 

2.8 

40039 

10 00 5 

-6022 

.760 8 

.0033 

.0076 

.6135 

D01S 

JD050 

2.9 

4*4 0 4 0 

10003 

-0014 

.761 0 

.0016 

.0055 

.6136 

D011 

-4)037 

3-0 

46040 

10003 

-0009 

.7611 

.0018 

.0040 

.613 7 

J0008 

4)027 

3.1 

4004 0 

ID 0 01 

-0 0 06 

.761 3 

.0008 

.0039 

.6138 

D 0 0 5 

D 0 19 

3.3 

4.7040 

ID 001 

-0 0 04 

.7613 

.0006 

.0031 

.6138 

D 0 0 4 

D014 

33 

40 04 0 

ID 0 01 

-0 003 

.7613 

.0004 

.0015 

.6139 

D 0 0 3 

D010 

3^4 

AS 0 4 0 

ID 0 0 0 

-D001 

.761 4 

.0003 

.0010 

.6139 

D 0 0 2 

4)007 

3.5 

5D040 

ID 0 0 0 

-oooi 

.761 4 

.0002 

.0007 

.613 9 

J0001 

.0005 

3.6 

53 0 4 0 

ID 0 0 0 

-oooi 

.761 4 

.00 01 

.0005 

•613 9 

.0001 

.0003 

37 

50040 

12000 

.0 0 00 

.761 4 

.00 01 

.0003 

*613 9 

.0001 

.000 2 

38 



.7614 

.0001 

.0002 

•6139 

D 0 0 0 

.0 0 02 

39 




.7614 

.0000 

.0002 

.6139 

JOOOO 

.00 01 

40 




.761 4 

.0000 

.0001 

.6139 

JO 000 

.0 0 01 

4J. 




.761 5 

.0000 

.0001 

.6139 

.0000 

.00 00 

4.3 

40 

4/# 




.761 5 

.0000 

.0000 

.6139 

J) ooo 

.0000 
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TABLE I. - Concluded. LAMINAR BOUNDARY LAYER BEHIND SHOCK WAVE (EQS. 9 AND 15) 
(g) Solution for u^Ug = 6 (ref. 1) 


17 

f 

f 


I Prandtl numbe? ow# 0.72 

r’l 

/ r dlj 

Jo 

r 

-r 1 

J 

8 

-8 « 

010 

04)0 0 0 

6D00 0 

-84. 0 09 

04)0 0 0 

0.9 195 

0.0 0 0 0 

0.00 0 0 

ID 0 00 

U 2 6 3 

.1 

,5597 

54977 

-72 7 31 

D 9 1 3 

.9 0 09 

•3658 

.0 934 

■8683 

1.2 991 

.a 

ID 410 

4 A 3 8 3 

-7D 6 2 0 

.17 91 

•8486 

•6668 

.1738 

.7417 

ID 2 5 8 

-3 

3,4 4 9 9 

3.7 53 4 

-6.4 078 

D6 02 

.7 710 

.8 671 

.24 20 

•6242 

1.12 0 3 

•4 

1.7 9 4 7 

34 6 0 3 

-5b4 4 5 4 

.3328 

.6 7 87 

.9 61 9 

.3991 

•5183 

S 963 

•5 

2D 8 5 2 

35643 

-4.4 8 31 

.3958 

.5 816 

D 6 8 1 

.34 61 

•4252 

D 663 

•6 

23 307 

2D 613 

-3.5937 

.449 2 

.4871 

.9118 

.38 4 5 

.3 4 50 

.73 8 7 

.7 

2J5 402 

ID 418 

-38161 

.4 93 5 

.4 0 04 

.8 187 

.4155 

.2771 

.6197 

•8 

2.7 2 1 4 

1.6 9 3 8 

-2.16 43 

.529 6 

.3239 

.70 9 4 

.44 0 3 

.32 06 

.5127 

.9 

2D 8 0 9 

15048 

-16353 

.558 6 

.2 5 86 

.5 98 3 

.4600 

•17 4 2 

.4190 

1.0 

3D 24 0 

13631 

-1.2171 

5817 

.2 04 0 

.4 941 

.4755 

.1364 

.3388 

1.1 

3.1 5 4 6 

ID 5 8 3 

-.8 9 3 5 

.599 8 

.15 94 

.4012 

.4 8 75 

.1060 

.3713 

1.2 

3D 7 6 6 

1 J. 8 I 8 

-6478 

.613 8 

.1234 

.3 212 

.4 9 69 

D 818 

.8151 

1.3 

33918 

14 26 6 

-.4 6 41 

.6 34 7 

.0 947 

j3 5 4 2 

.50 4 0 

D 6 26 

4.693 

3U4 

3D 02 4 

ID 8 7 3 

-.3 2 8 8 

.633 0 

.0 731 

4.991 

.50 9 5 

4)476 

.13 3 0 

1.5 

36097 

1.0 5 97 

-.3 3 0 4 

.6 39 3 

.0 545 

.15 4 5 

.5137 

.0 360 

.10 3 3 

1.6 

3.7 1 4 6 

ID 40 3 

-.159 7 

.6 44 0 

■0 409 

1189 

.5168 

.0 270 

.0785 

1.7 

3D 18 0 

ID 370 

-.109 6 

.6 47 6 

.0 3 05 

.0 908 

.5192 

.0 201 

.0 59 9 

3-8 

3D 20 3 

ID 179 

-.0 7 44 

.6 50 2 

.0 22 6 

.0 6 87 

.53 0 9 

D 1 4 9 

.0 4 53 

1-9 

44)317 

ID 11 8 

-.0 5 00 

.6521 

.016 6 

.0517 

.52 2 2 

D 109 

.0340 

2.0 

44 32 6 

ID 077 

-.0 333 

.653 6 

.0131 

.0385 

.5331 

.0 0 80 

.035 4 

3.1 

4.2 333 

ID 04 9 

-D 219 

.654 6 

.0088 

.0 285 

.52 3 8 

D 0 5 6 

.0188 

2.2 

43336 

14)031 

-D 14 3 

.655 4 

.0 06 3 

.0210 

.53 4 3 

.0042 

.0138 

2.3 

4.4 239 

1-0 0 3 0 

-D 09 2 

.655 9 

.0 04 5 

.0153 

.52 4 6 

.0 0 30 

.0101 

2.4 

4Ji 34 0 

ID 01 2 

-D 0 5 9 

.656 3 

.0 032 

.0111 

.5249 

.0 0 31 

.00 7 3 

2^ 

4.6 341 

ID 0 08 

-D 0 37 

.656 5 

.0033 

.00 80 

.53 51 

.0015 

.00 5 3 

2.6 

4.7 2 4 2 

14)005 

-DO 23 

.656 7 

.0 016 

.0 0 57 

.52 5 2 

.0 010 

.00 3 8 

a.7 

4.8 24 2 

10 00 3 

-D 015 

.656 9 

.0011 

.0 04 0 

.5253 

.0 0 07 

.00 3 7 

2.8 

49343 

1.0 0 0 2 

-.0 0 09 

.6570 

.0 0 08 

.00 28 

.5253 

•0 0 05 

.00 19 

2.9 

SL0343 

10001 

-DO 05 

.657 0 

.00 05 

.0020 

.5254 

•0003 

.0013 

3.0 

31343 

1.0001 

-DO 03 

.6571 

.0004 

.0014 

.5354 

•0002 

.00 0 9 

3.1 

5J3 34 3 

10000 

-0 0 03 . 

-6571 

.0003 

.0010 

.5354 

.0 0 02 

.0006 

3.3 

33 84 3 

10000 

-.0001 

.6571 

.0 0 03 

.0007 

.5354 

.0001 

.0 00 4 

3.3 

34343 

10000 

-4) 0 01 

.6571 

.0001 

.000 4 

.5255 

.0 0 01 

.0003 

3.4 

35 34 3 

10000 

.0000 

-6571 

.0001 

.0003 

.5355 

.0 0 00 

.0003 

35 




.6572 

.0000 

.0002 

.5255 

.0 0 00 

.0001 

3.6 




.657 3 

.0000 

.oooi 

.5355 

.0 0 00 

.0 0 01 

37 




.6573 

.oooo 

.0 0 01 

.5355 

.0 0 00 

.0001 

38 




.6573 

.0 0 00 

.0001 

.5255 

.0 0 00 

.00 0 0 

39 




.6572 

.0000 

.00 00 

.5255 

.0 0 00 

.0000 


3959 









3959 
























TABLE III. - COMPARISON OF INTERPOLATION FORMULAS WITH EXACT INTEGRATIONS 
OF EQUATIONS (9) AND (15) (FRANDTL HUMBER <%, 0.72) 




0.75 




r r dT) 

J o 


Tad, 

^0 


- 8 ’( 0 ) 


r(0) 


Eqs. 

(9), (15) 

Eg.. 

(27) 

Eqs. 

(9), (15) 

Eg. 

(27) 

Eqs. 

(9), (15) 

Eg. 

(27) 

Egs. 

(9), (15) 

Eg. 

(27) 

Egs. 

(9), (15) 

Eg. 

(27) 

0.5727 

0.570 

0.6574 

0.655 

0.7312 

0.730 

1.0191 

1.017 

1.3541 

1.353 

1.0448 

1.054 

0.9360 

0.944 

0.8580 

0.862 

0.6466 

0.650 

0.5013 

0.503 

1.3970 

1.396 

1.2822 

1.279 

1.1964 

1.194 

0.9482 

0.940 

0.7615 

0.754 

1.1974 

1.201 

1.0854 

1.087 

1.0038 

1.004 

0.7767 

0.774 

0.6139 

0.610 

0.4996 

0.502 

0.5665 

0.569 

0.6247 

0.627 

0.8512 

0.856 

1.1156 

1.122 

0.8628 

0.858 

0.8727 

0.867 

0.8799 

0.876 

0.8997 

0.893 

0.9129 

0.907 
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TABLE 17. - RECIPROCALS OF OTEEGRALS OCCURRING IN TURBULMT-BOORBARr-LAZER ANALYSIS 
(a) Reciprocal of Ig (based on numerical results of ref. ll) 



1.958 

1.830 

1.700 

1.567 

1.430 

1.288 

1.139 

.9814 

.8969 

.8071 

.7097 

. 59 % 

.4575 

.3469 

0 


3.270 

3.036 

2.799 

2.557 

2.310 

2.056 

1.792 

1.516 

1.370 

1.217 

1.054 

.8726 

.6418 

.5003 

0 


-.4 

-.2 

0 

4.531 

5.770 

6.998 

4.193 

5.328 

6.452 

3.851 

4.883 

5.900 

3.503 

4.427 

5.342 

3.149 

3.967 

4.776 

2.784 

3.495 

4.196 

2.410 

3.011 

3.604 

2.019 

2.507 

2.988 

1.814 

2.244 

2.716 

1.599 

1.970 

2.334 

1.372 

1.682 

1.983 

1.123 

1.366 

1.603 

.8285 

.9960 

1.161 

.6317 

.7559 

.8764 

0 

0 

0 


19.07 

17.49 

15.92 

14.31 

12.70 

11.06 

9.390 

7.680 

6.789 

5.869 

4.921 

3.897 

2.738 

2.027 

0 


31.00 

28.43 

25.83 
23.19 
20.54 

17.84 
15.11 
12.29 
10.86 

9.363 

7.831 

6.150 
4.288 

3.151 
0 


42.96 

39.34 

35.74 

32.06 

28.39 

24.63 

20.79 

16.91 

14.89 

12.84 

10.70 

8.410 

5.865 

4.270 

0 


54.88 

50.13 

45.62 

40.90 
36.17 
31.43 

26.52 

21.52 

18.91 
16.31 
13.58 

10.63 
7.391 
5.397 
0 


(b) Reciprocal of 1^ (based on numerical results of ref. ll) 



2.175 

2.027 

1.879 

1.726 

1.570 

1.409 

1.240 

•1.062 

.9672 

.8668 

.7584 

.6360 

.4805 

.3606 

0 


3.683 

3.411 

3.135 

2.855 

2.570 

2.278 

1.976 

1.660 

1.495 

1.522 

1.138 

.9350 

.6924 

.5247 

0 


-.4 

-.2 

0 

5.142 

6.581 

8.009 

4.746 

6.061 

7.366 

4.345 

5.536 

6.717 

3.940 

5.006 

6.062 

3.528 

4.46B 

5.398 

3.107 

3.920 

4.724 

2.675 

3.359 

4.034 

2.226 

2.778 

3.322 

1.992 

2.476 

2.952 

1.748 

2.163 

2.571 

1.491 

1.834 

2.170 

1.211 

1.478 

1.739 

.8826 

1.064 

1.243 

.6657 

.7987 

.9276 

0 

0 

0 


22.12 

20.25 

18.36 

16.47 

14.56 

12.62 

10.66 

8.643 

7.609 

6.548 

5.447 

4.278 

2.965 

2.165 

0 


36.16 

33.06 

29.94 

26.81 

23.65 

20.46 

17.22 
13.91 

12.22 

10.48 

8.685 

6.783 

4.662 

5.381 

0 


50.18 

45.85 

41.50 

37.13 

52.73 

28.28 

23.78 

19.17 

16.82 

14.41 

11.91 

9.282 

6.354 

4.593 

0 


64.20 

58.64 

53.06 

47.45 

41.80 

36.10 

30.32 
24.43 
21.41 

18.33 
15.14 
11.78 

8.044 

5.804 

0 


(c) Reciprocal of Ig (based on numerical results of ref. ll) 



2.388 

2.221 

2.054 

1.882 

1.707 

1.527 

1.338 

1.140 

1.035 

.9242 

.8045 

.6711 

.5023 

.3734 

0 


-.6 

-.4 

-.2 

0 

4.092 

5.744 

7.380 

9.001 

3.781 

5.294 

6.784 

8.264 

3.467 

4.836 

6.184 

7.519 

3.149 

4.373 

5.574 

6.770 

2.827 

3.905 

4.963 

6.010 

2.496 

3.425 

4.558 

5.241 

2.156 

2.936 

3.701 

4.456 

1.801 

2.430 

3.045 

5.650 

1.617 

2.166 

2.703 

3.231 

1.424 

1.893 

2.351 

2.801 

1.219 

1.606 

1.983 

2.352 

.9950 

1.295 

1.586 

1.871 

.7294 

.9337 

1.129 

1.322 

.5473 

.6974 

.8389 

.9766 

0 

0 

0 

0 


25.11 

22.95 

20.76 

18.58 

16.38 

14.15 

11.90 

9.597 

8.418 

7.210 

5.963 

4.645 

3.179 

2.296 

0 


41.17 

37.57 

33.96 

30.33 

26.69 

23.00 

19.28 

15.48 

13.56 

11.57 
9.542 
7.386 
5.010 
3.592 
0 


57.18 

52.16 

47.10 

42.07 

36.98 

31.85 

26.67 
21.38 

18.68 
15.93 
13.09 
10.12 

6.840 

4.885 

0 


73.21 

66.76 

60.28 

53.82 

47.26 
40.68 
34.03 

27.26 
23.80 
20.28 
16.66 
12.86 

8.666 

6.177 

0 













Figure 2. - Temperature and velocity distributions behind a shock wave. 


